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Abstract 

Large back-to-back correlations of observable fermion - anti-fermion pairs are 
predicted to appear, if the mass of the fermions is modified in a thermalized 
medium. The back-to-back correlations of protons and anti-protons are exper- 
imentally observable in ultra-relativistic heavy ion collisions, similarly to the 
Andreev reflection of electrons off the boundary of a superconductor. While 
quantum statistics suppresses the probability of observing pairs of fermions 
with nearby momenta, the fermionic back-to-back correlations are positive 
and of similar strength to bosonic back-to-back correlations. 
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I. INTRODUCTION 



In high energy heavy ion collisions, hot and dense hadronic matter is expected to be 
created in conditions similar to the ones in the early Universe about a few /isec after the 
Big Bang. Under these conditions, strongly interacting particles may propagate with a mass 
that differs from the mass in the asymptotic vacuum. Recently, it has been discovered that 
in-medium mass-modifications induce large back-to-back correlations (BBC) among pairs of 
asymptotic, observable bosons However, it was not known before if the bosonic BBC 

would have an analogous effect in the fermion sector. In this Letter we show that fermion - 
anti-fermion pairs also exhibit BBC, if they propagate with a modified mass in the medium. 
The physical picture is that of a system in thermodynamical equilibrium which freezes-out 
suddenly. The masses and other properties of the hadrons in the system are modified due to 
interactions and a quasi-particle description is assumed. At freeze-out the quasi-particles are 
converted into the free particles that will be detected. Particle-anti-particle pairs detected 
with momenta back-to-back give rise to a quantum-mechanical correlation due to a nonzero 
overlap between the in-medium states and the free states. The quantum correlation can be 
calculated with a Bogoliubov-Valatin transformation between a quasi-particle basis and a 
free-particle basis. 

The bosonic BBC (bBBC) have a quantum optical analogy, namely the correlations in 
thermalized ensembles of two-mode squeezed states. It turns out that fermion BBC (fBBC) 
also have an analogy, which is the Andreev reflection, well-known in solid state physics. It 
refers to the scattering of electrons off the boundary of a superconductor - normal conductor 
junction. The reflected electrons are used to study the properties of the superconductors || . 
Our results for fBBC generalize Andreev's reflection to the case when superconductivity is 
suddenly switched off in the whole volume of the material, so that a junction prevails not 
at a given position for a long time but in the whole medium at a given instant. 

Quantum statistics enhances the probability of observing pairs of bosons in similar mo- 
mentum states, while it suppresses the probability of observing pairs of fermions with nearby 
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momenta. In spite of such Fermi-Dirac suppression factors that appear in the evaluation of 
the correlation function, we find that the dominant term of the fBBC, due to the Bogoliubov- 
Valatin transformation, is positive and similar in strength to the bosonic BBC of Refs. |I]|§. 
The fBBC could be observed experimentally in 197 Au + 197 Au collisions at the Relativistic 
Heavy Ion Collider (RHIC), which started to take data at the Brookhaven National Labo- 
ratory at y/s = 56 and 130 AGeV in 2000 and which will reach its full designed energy of 
y/s = 200 AGeV in the near future. 

The outline of the paper is as follows. In section II we consider the in-medium mass 
modification state at finite temperature. The spectra and the correlations for mass-shifted 
fermions are described in section III and section IV contains concluding remarks. 

II. BASIC ASSUMPTIONS 

We assume the validity of the concepts of local thermalization, hydrodynamics, and a 
short duration of particle emission. These concepts are in agreement with the observable 
single-particle spectra and two-particle correlations of pions, kaons and protons [|]]. We also 
assume the validity of the effective Hamiltonian 

H = H + Hj, (1) 

where 

H = J rfx : $(x)(-i7 ■ V + M)^(x) : (2) 

is the Hamiltonian in vacuum, and the interaction Hamiltonian Hj describes the medium 
modifications and will be discussed shortly below. In Eq. @, M is the value of the proton 
(or in general, baryon) mass in free space, ip and ip are the fermion field operators which 
satisfy the usual equal-time anti-commutation relations. 

We are interested in the single-particle and two-particle invariant momentum distribu- 
tions for fermions and anti-fermions: 
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iVi(ki) = ^(a^aki), iVi(ki) = a; kl (a kl a kl ), (3) 
iV 2 (ki,k 2 ) = ^^(al^a^a^). (4) 

Here, (O) denotes the expectation value of operator O in the thermalized medium, and the 
a', a, a', a are creation and annihilation operators of free baryons and anti-baryons of mass 



M and energy uj k = \J M 2 + |k| 2 . These creation and annihilation operators are defined 
through the expansion of the baryon field operator as 

^( x ) = —fff ( u A,kaA,k + v\> _k4' ,-k) e * k ' X ' ( 5 ) 
W \,\',k 

where V is the volume of the system, and t>A',- k are Dirac spinors and the summation 
extends over momenta k and spin projections A, A' = 1/2,-1/2. The creation and anni- 
hilation operators satisfy canonical anti-commutation relations. The particle-anti-particle 
correlation function is defined as 

C 2 (k 1 ,k 2 )= N2{klM) . (6) 
21 V iV 1 (k 1 )iV 1 (k 2 ) 1 ; 

As in the bosonic case , the expectation value of the four baryon operators in Eq. (|J) 
can be calculated using Wick's theorem generalized for a locally equilibrated system, with 
the result 

(akA^aki) = (4i a ki)(Gk 2 5 k 2 ) - (a ki ak 2 )(4 2 a ki) + ( a k 1 4 2 )(«k 2 ak 1 )- (7) 

The minus sign in the above equation is due to Fermi statistics. The expectation values 
involving a ki ak 2 and a kl 5k 2 give rise to the chaotic amplitudes, while a kl a k2 gives rise to 
the "squeezed" amplitude, defined as 

G c (ki,k 2 ) = v /^kT^k 2 "(a kl «k 2 ) , G c (ki,k 2 ) = A /w k ^ k 7(a kl o k2 ), (8) 
G s (ki,k 2 ) = ^/uJ^O^{a kl a k2 ) . (9) 

The chaotic amplitudes are non-vanishing for identical type of fermions and anti-fermions, 
while the squeezed amplitude is non-vanishing only for particle-anti-particle pairs. We will 



show in the following that, as in the case of bosons [0,0, the squeezed G s (ki,k 2 ) can be 
considerably large if the in-medium masses of the baryons are different from their free-space 
values. 

In order to evaluate the thermal averages above we model the system as a globally 
thermalized gas of quasi-particles (quasi-baryons). The density matrix for such a system 
is given not in terms of the free-space creation and annihilation operators a', a, ■ ■ ■ , but 
by a new set that we denote by , b, • • •. They are of course different because of medium 
effects, described by the interacting Hamiltonian in Eq. ([!]). While it is the a-quanta that 
are observed as asymptotic states, it is the b quanta that are thermalized in the medium. 

In a quasi-particle description of the system, the medium effects are taken into account 
through a self-energy function. For a spin-1/2 particle under the influence of mean fields in 
a many-body system, one can write its self-energy function as 

E = S s + 7°S° + 7*S i . (10) 

This function can be determined from a detailed calculation based on, for example, the 
Walecka or the Zimanyi - Moszkowski models |S]||. Calculations in the context of such 
models show that E° is weakly momentum dependent and is very small (see for example 
Ref. H for detailed numerical results). It is well known that a tensor contribution may also 
occur, however, in applications for nuclear matter and finite nuclei this additional component 
is negligible, and will not be considered here. Therefore, we assume 2° independent of 
momentum, and neglect E\ We denote the scalar component of the self-energy as S s (k) = 
AM(k). Now, for a locally thermalized system, the role of S° is to shift the chemical 
potential, i.e. 

At* = M-£°. (11) 

As we present our results as functions of the net baryon density, the value of S° (or the 
difference between \i and /z*) and its \x dependence need not be specified. More specifically, 
since we specify the total baryon density, we invert the expression of the density in terms of /x* 
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to obtain the baryon and anti-baryon Fermi-Dirac distributions (see Eq. fll7|) below). In this 
way, the value of S° and its dependence on \i are not needed. With these approximations, the 
effective Hamiltonian of Eq. ([[]) thus describes a system of quasi-particles with a momentum- 
dependent mass M*(k) — M — AM(k). 

In the context of the model we just described, the free-space and in-medium creation and 
annihilation operators are related through a fermionic Bogoliubov-Valatin transformation @ 

\ 



/ a \ ( 

a A,k 



a A',-k 



Ck 



fos k A\( b - ^ 



l/i 



V i/k 



/ 



>»A,k 



St 



A' — k 



(12) 



where Ck = cos/k, Sk = sin/ k , ^4 is a 2 x 2 matrix with elements A\ t \i = Xx 17 ' ^Xx'f where 
k = k/|k|, x is a Pauli spinor and \ = —io 2 \. Here, /k is the squeezing function. In the 
present case, f k is real and we will therefore drop the complex-conjugate notation in what 
follows. The effective Hamiltonian is brought to a diagonal form in the basis of the b quanta 
when the squeezing function f k is related to the mass shift AM as [§], 

|k|AM(k) 



tan(2/ k 



(13) 



ojI - MAM(k) ' 

With this, the thermal averages are then evaluated with the following density matrix oper- 
ator 



1 



P 



exp 



1 V 



dk 



(14) 



Z ~ r V T (2tt) 3 

where V is the volume of the system, T the temperature, and f2k = (k 2 + M* 2 ) 1 / 2 is the 
quasi-particle energy. Z is the trace of the exponential factor in Eq. ( |T1D . 

The evaluation of the thermal averages in Eqs. dg) and (^) is very easy and the results 
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( a A,k fl A',k') = (4 n k + 4 (1 - n k )) S X x' S kk >, 
(a x - k a X 'y, ) = (1 - n k - n k )c k s k (A)aA' <W 



(15) 
(16) 



where 
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exp [(fi k - ^)/T] + 1 ' exp [(O k + /i*)/T] + 1 ' y J 

The terms (a Ak a,y,k') and (a^ k aA',k') in Eq. (^) vanish because they involve expectation 
values of two baryon annihilation operators, or two anti-baryon creation operators, or one 
baryon and one anti-baryon operators. 



III. SPECTRA AND CORRELATIONS FOR MASS-SHIFTED FERMIONS 



If a thermal gas of b fermions freezes out suddenly at temperature T, the observed single 
particle and anti-particle spectrum are given by 

V 



iVi(k) 



jVi(k) 



(2k)- 



V 



c£n k + 4(1 - n-k) 



4(1 - n k ) + c k n k 



18a) 



;i8b) 



Note that the single-particle (anti-particle) spectrum includes a squeezed contribution. 
While the thermal part of the spectrum falls off exponentially for large values of |k|, the 
squeezed contribution falls off only as a power of |k|, because of the term proportional to s k 
in Eq. fllf). 

For such a homogeneous system, the net baryon density p#, the chaotic, and the squeezed 
amplitudes can be written from Eqs. ( ]T5|) and ([16]). They are given by 



PB = TfZI (™k - *»k), 

v k 



G c (k!,k 2 ) 
G c (k ls k 2 ) 



V 



(2tt) 3 
V 

V 
J2nf 



Wki [4i "ki + 4i (! - »*ki)] <W 2 , 
^kx [4i (1 - "ki) + 4x ™k x ] <W 2 , 
^ki [(1 - n kl - n kl )c kl s kl A f ] 5 kl ,-k 2 • 



(19) 
(20) 
(21) 
(22) 



There are two kinds of fermionic two-particle correlation functions. Similarly to the bosonic 
case 0, let us denote by C^" 1 " (ki, k 2 ) the case when the two particles are identical fermions, 
and by C 2 + ^(ki, k 2 ) when particle 1 is a fermion and 2 is an anti-fermion. The correlation 



functions C2 (kx, k 2 ) and C2 (k 1; k 2 ) can be obtained from the (++) and (H — ) correla- 
tions by a trivial exchange of particle and anti-particle labels. For an infinite, homogeneous 
thermalized medium, these correlation functions are non-trivial only for identical or back-to- 
back momenta, k 2 = ±ki. We find that C*2 ++ ' ) (k, k) = C% (k, k) = 0, the canonical value 
of the Fermi-Dirac correlation function, reflecting the anti-correlation of identical fermions, 
due to the Pauli exclusion principle. 

The BBC for fermion - anti-fermion pairs given by Eq. (^), reads as 

|C,(k, -k)l 2 



Cf- } (k,-k) = 1 + 



G c (k,k)G c (-k,-k) 



= x + (1 -n k -n k ) 2 (c k g k ) 2 

[cW + 4(1 - n k )] [c k n k + 4(1 - n k )] ' 

From this equation it follows that the fermionic BBC are unlimited from above. For suf- 
ficiently large values of |k|, the Fermi-Dirac distribution n k falls exponentially, while s k 
decreases only as a power- law. Hence, for sufficiently large values of |k|, the fBBC diverge 
as in the case of bosonic BBC, as \k, — k) oc 1 + l/s k — > 00. This divergence happens 
for small values of the mass-shift and for large values of k in both the fermionic and the 
bosonic cases. 

We next discuss numerical results. We use momentum-independent in-medium masses, as 
one would obtain in a mean-field approximation in model calculations of the sort discussed 
in Refs. |5]||. There is no intrinsic difficulty in using momentum-dependent self-energies 
but a commitment to a specific model is then necessary. Since in this initial study we are 
more interested in the qualitative effects, rather than in precise, quantitative predictions, we 
simply use typical values for M* as obtained in mean-field calculations. In Fig. 1 we show the 
fBBC for pp pairs as a function of the in-medium mass M* = M — AM, for three illustrative 
values of the net baryonic density p#: the normal nuclear matter density, one tenth of the 
nuclear matter density, and p% = 0. This last value of pB corresponds to a baryon-free 
region, as expected to be formed at RHIC. We also show in this figure the corresponding 
result for bBBC of mesons. Following Ref. 0], a finite time suppression factor is used to 
model a more gradual freeze-out. A detailed derivation for a gradual exponential freeze-out 
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leads to a factor of the form , 

^ (k > |2 = [i + (2L k )T (24) 

which multiplies C2 (k, k) — 1, the non-trivial part of the back-to-back correlation function. 
We use At = 2 fm/s, as in Ref. ||. We observe that the fBBC are strongly enhanced when 
the net baryonic density decreases, and that the shape of fBBC for ps = is rather similar 
to bBBC. For a fixed temperature, as we increase the net baryon density obviously we 
have an excess of baryon over anti-baryons. On the other hand, the fBBC will be larger 
for approximately equal baryon and anti-baryon densities and therefore the effect will be 
enhanced as we approach zero net baryon density. Our result shows that for a baryon free 
region, the BBCs are approximately independent of the bosonic or the fermionic nature of 
the particles. In addition, the fBBC are not only positive, as in the bosonic case, but they 
are also of the same order of magnitude as the corresponding bBBC. 

The expected magnitude of the effect is illustrated in Fig. 2 for two, typical momenta 
of thermal-looking proton spectrum in Pb + Pb collisions at CERN SPS 0, for |k| = 500 
MeV and for |k| = 800 MeV. As in Fig. 1, the finite time suppression factor of Eq. ( |24D 
is used. Fig. 2 indicates that the fBBC is strongly enhanced for increasing momentum |k|. 
Fig. 2 also highlights that the magnitude of the fBBC is greatly enhanced as the net baryon 
density decreases from normal nuclear density to a vanishing value. Similarly to the bBBC 
case, the strength of fBBC is very sensitive to the shape of the freeze-out distribution, as 
well as to the value of the parameter At, the particle freeze-out duration, as discussed in 
Ref. 0. Note that in the limit of infinitely slow freeze-out, At — > oo, both fBBC and bBBC 
vanish. 



It is particularly interesting to compare the BBC discussed above in Eq. (|23|) with 
Eq. (19) of Ref. Q for bosonic 

cf~\K -k) = l + lt ^.. t ^\^}±^\^. ^ , (*) 



(4)2 ni + fe) 2 (1 + ni) (4) 2 ni + (4) 2 (1 + v h 
1 

exp (Q b jT) 



< = ^ZTTTWl^ i • ( 26 ) 



Note, however, that the symbols Sk and Ck are related to sin/k and cos/k, as defined by 
Eq.(|l3D in the fermionic case, while in the bosonic case, these symbols mean s k = sinh/^ 
and c kB = cosh/k, where /£ = \ log (c^k/^k), as given by Eq. (11) of Ref. [TJ], is also real. It 
follows that the BBC diverge with increasing |k| as the inverse of the single-particle spectra, 
both in the fermionic and bosonic cases. 

IV. CONCLUSION 

BBC stands for Back-to-Back Correlations of particle and anti-particle pairs. These 
correlations appear if in-medium interactions lead to the modification of hadronic masses 
in the medium. The origin of these back-to-back correlations is an entirely quantum effect, 
related to the propagation of particle fields through a space-like boundary surface between 
the medium and the asymptotic region. The strength of these correlations can be unlimitedly 
large, and the shape of the BBC is similar for fermions and bosons. Deep mathematical and 
physical reasons are behind these similarities. A sudden freeze-out of thermalized medium- 
modified quanta to asymptotic fields is described by a bosonic or a fermionic Bogoliubov 
transformation. Although these transformations are canonical, they connect Fock spaces 
that become unitarily inequivalent in the infinite volume limit and the state of the medium 
corresponds to strongly correlated, squeezed particle-anti-particle states of the asymptotic 
quanta. 

At large values of momenta, the similarity of fermionic BBC to bosonic BBC reflects 
a symmetry of the decaying medium to observable boson - anti-boson and fermion - anti- 
fermion pairs. This effect vanishes in case of exactly zero in-medium mass-modification, but 
for large values of |k| and small values of mass-shifts, BBC can be very large and should be 
observable if the freeze-out process of the medium modified quanta is sufficiently sudden, 
for example if the duration of particle freeze-out is within 2 fm/c. 

In this Letter, we have extended the concept of back-to-back correlations to an important 
new domain of broad experimental and theoretical interest. As a reduction of the net baryon 
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density by a factor of 10 increases dramatically the magnitude of the effect, the almost 
baryon-free mid-rapidity region of Au + Au collisions at RHIC seems to be an ideal place to 
find the back-to-back correlations of proton - anti-proton or A - A pairs experimentally. 
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FIG. 1. Back-to-back correlations of proton - anti-proton pairs and <^>-meson pairs, for 
T = 140 MeV, At = 2 fm/c and |k| = 800 MeV/c. 
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k=800 MeV 




FIG. 2. Dependence of the fBBC on the in-medium modified proton mass, M*, and on the net 
baryon density, ps, for T = 140 MeV, At = 2 fm/c and two typical values of the momentum |k|. 
The fBBC strongly increases with decreasing net baryon density and increasing values of |k|. 
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